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1 Introduction 



It is well established that density inhomogeneities play an important role in 
the collapse of a spherical dust cloud. In particular they may lead to the 
formation of naked singularities H, in contrast with the homogeneous case 
leading to a black hole 0. However, it is unclear what is the physical reason 
(if any) behind the link between the final fate of collapse and the presence 
(or absence) of density inhomogeneities. 

It is the purpose of this work to explore the ways in which density inho- 
mogeneities may affect the spherical collapse. To do so, we shall obtain an 
expression for the active gravitational mass (the Tolman mass) containing 
explicitly a measure of density inhomogeneity (among other factors). Since 
our expression is obtained for a general (locally anisotropic) fluid, it also 
contains a term depending on local anisotropy. We shall see that this term 
plays a similar role to the density inhomogeneity term, thereby suggesting 
that local anisotropy might also lead to the formation of naked singularities. 
As a by-product of this study we shall find some expressions linking the Weyl 
tensor and the mass function to the energy density inhomogeneity and local 
anisotropy. 

The paper is organized as follows. In the next section the field equations and 
other useful formulae are introduced. In section 3 we introduce the Tolman 
mass and discuss its physical meaning. In Section 4 we derive an expression 
for the Tolman mass and evaluate it at the moment the system departs from 
hydrostatic equilibrium. A discussion of this expression is presented in the 
last section. 

2 Field Equations and Conventions 

We consider a spherically symmetric distribution of collapsing fluid, which 
for completeness we assume to be anisotropic and bounded by a spherical 
surface S. The line element is given in Schwarzschild-like coordinates by 

ds 2 = e v dt 2 - e x dr 2 - r 2 (dO 2 + sin 2 6d<f) 2 ) (1) 

where v and A are functions of t and r . The coordinates are: x° = t; x 1 = 
r; x 2 = 8; x 3 = 0. 

The metric ([I]) has to satisfy Einstein field equations 
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G$ = -MT2 (2) 

which in our case read HI: 
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-8ttT 01 = -- (6) 
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where dots and primes stand for partial differentiation with respect to t and 
r respectively. 

In order to give physical significance to the components we apply the 
Bondi approach 0, i.e we introduce local Minkowski coordinates (r,x,y,z), 
defined by 



dr = e u l 2 dt dx = e x ^ 2 dr dy = rd6 dz = rsin 

Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



T ° = T ° Tl = T\ Tl = T| T 3 3 = T 3 T 01 = e - {u+x)/2 T m 

Next we suppose that, when viewed by an observer moving relative to these 
coordinates with velocity uj in the radial direction, the physical content of 
space consists of an anisotropic fluid of energy density p, radial pressure P r 
and tangential pressure P±. Thus, when viewed by this moving observer, the 
covariant energy-momentum tensor in Minkowski coordinates is 
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Then a Lorentz transformation readily shows that 

- p + P r cu 2 
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Note that the velocity in the (t, r, 6, 0) system, dr/dt, is related to by 

w = d J_ e (A-,)/2 

Outside of the fluid , the spacetime is Schwarzschild, 



(7) 

(8) 
(9) 

(10) 
(11) 



- - ' 1 dl- 
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ds 2 = 1 



1 _ ^ 1 rfr 2 _ r 2 ^2 + sm 2^ 2^ (12) 



In order to match the two metrics smoothly on the boundary surface r = 
rs(i), we require continuity of the first and second fundamental forms across 
that surface. As result of this matching we obtain the well known result 







The radial component of the conservation law 



= 



(13) 
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gives 
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(-8<)' = ^ (7? _ jj) + W (7? - T°) + ^ f A + f - ^) (15) 



which in the static case becomes 



P , r = ~(p + Pr)+ 2{Pl ; Pr) (16) 

representing the generalization of the Tolman-Oppenheimer-Volkof equation 
for anisotropic fluids 

For the next sections it will be useful to calculate the components of the 
Weyl tensor. Using Maple V, it is found that all non-vanishing components 
are proportional to 



W = W { . ) + ^(\ + £-%) (17) 



where 



Ww = — {^~^ + — ~T-Y~Yr + Yr) (18) 

corresponds to the contribution in the static case. 

Finally, let us introduce the mass function, defined by || |6| 

m(r } t) = l -rR^ 2 (19) 
where the Riemann component for metric ([I]) is given by 

^232 = l~e- x (20) 
Then using the field equation (||) we obtain the well known expression 

m(r,t) = 4vr F r 2 T °dr (21) 

J 

Now, considering the definition of the Weyl tensor 
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(22) 

and eqs. ©, §), § 
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171) and ([19]), it follows that 
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r 3 (T ° + Tl - T%) + W 



Differentiating (^) with respect to r and using ( 21 ) we get 



An 



An 



and integrating 



Finally, inserting (|25| ) into (p^ ) we obtain 

An ,„ n 47T 



4tt 



W = / r 3 (T °) dr + — r 3 (T| - 2? 



m r. 



t) = fr 3 T °-|/ o V 3 (r o ) , rfr 



(23) 



(24) 



(25) 



(26) 



Observe that expressions (|23"D, (|25| ) and fl26|) are the same for the static (or 
slowly evolving) case 0. Of course, in the latter case Tq and — Tf denote 
the proper energy density and radial pressure, which is no longer true in the 
general case as can be seen from eqs.(^) and (|]). 

Instead of dealing with the mass function we shall now use the Tolman mass, 
which, as will be seen in the next section, appears to embody the physical 
idea of active gravitational mass better than the mass function. 



3 The Tolman mass 

The Tolman mass for a spherically symmetric distribution of matter is given 
by (eq.(24) in §) 

m T = An jf* r 2 e^ 2 (t o ° - 7* - 2T|) dr 

2io dt[d[d(!f*y/=g)/dt]) 9 1 ' 
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where L denotes the usual gravitational lagrangian density (eq.(10) in H). 
Although Tolman's formula was introduced as a measure of the total energy 
of the system, with no commitment to its localization, we shall define the 
mass within a sphere of radius r, completely inside S, as 



m T = 4vr f r 2 e^ +x ^ 2 (t o ° - 7* - 2T 2 2 ) dr 

+ I /"V C +A)/2 — f — ^ a a/3 dr (28) 

This extension of the global concept of energy to a local level || is sug- 
gested by the conspicuous role played by vht as the "effective gravitational 
mass", which will be exhibited below. Even though Tolman's definition is 
not without its problems || [TDH , we shall see that my, as defined by (|28f) , 
is a good measure of the active gravitational mass, at least for the system 
under consideration. 

Let us now evaluate expression (^). The first integral in that expression 

I = 4tt f r 2 e ( » +x)/2 (t o ° - T\ - 2T 2 ) dr (29) 

may be transformed as follows. Integrating by parts and using fl2~I|) , we 
obtain 



>+A)/2 



47T 



m(r, t) - —r 3 [T^ + 2T, 
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dr + 



+ 
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T 1 1 )' + 2(T 2 2 )' 



dr 



(30) 



Using the field equations and (|15|) we then obtain 



>+A)/2 



,(A-iO/2^ 




(31) 
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From (eq.(13) in §) 



dt(d[d (g^^9) /dt] ) " T ° af3+ 2 6 ° aT ^ + 2^ L (32) 
and so the second integral (II) in (|28f) may be expressed as 

"=^( s+ t-t)* (33) 

Thus 

m T = I + II = e {u+x)/2 [m(r, t) - 4 7 rr 3 T 1 1 ] (34) 

This is, formally, the same expression for my in terms of m and that 
appears in the static (or quasi-static) case (eq.(25) in [[/J). 
Replacing by @ and m by (|T9|) and (p0|) , one also finds 

mT = e ^m u ' r l (35) 

This last equation brings out the physical meaning of as the active grav- 
itational mass. Indeed, it can be easily shown [llj that the gravitational 



acceleration (a) of a test particle, instantaneously at rest in a static gravi- 
tational field, as measured with standard rods and coordinate clock is given 
by 

e (u-X)/2 I 

(36) 



2 r 2 

A similar conclusion can be obtained by inspection of eq . (|i~6|) (valid only in 
the static or quasi-static case) ]I2|. In fact, the first term on the right side of 
this equation (the "gravitational force" term) is a product of the "passive" 
gravitational mass density (p + P r ) and a term proportional to m-r/r 2 . 
In the next section we shall get another expression for m^, which appears 
to be more suitable for the treatment of the problem under consideration. 
This expression will be evaluated immediately after the system departs from 
equilibrium, where "immediately" means on a time-scale such that the value 
of uj remains unchanged. Therefore the physical meaning of m-r, as the 
active gravitational mass obtained for the static (and quasi-static) case, may 
be safely extrapolated to the non-static case within this time-scale. 
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4 Tolman mass, density inhomogeneity and 
local anisotropy 

The required expression for the Tolman mass will be obtained as follows. 
Taking the r-derivative of (^) we obtain 



>-A)/2 



rm n 



v' 2 AV „ 2v' 

h v H 
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On the other hand, it follows from (153) and fl23T) that 



(37) 



3m T = e [u+x) [4vrr 3 (t o ° - 27? - T|) + 3W] (38) 
Combining (|3?|) and (p8|) and using the field equations and (113), we obtain 



rm'rp — 3rriT = e 
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(39) 



which can be formally integrated to give 
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or equivalently 
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Finally, evaluating (|0) at to obtain C, and replacing it in (|T), we obtain, 
using (0) and (^5|) , 



m T = (m T )s (J^J 

- r 3 r s e ( ^ +A)/2 



5 7T 



r 



Z? - T. 



r 4 Jo 
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(42) 



In the static (or quasi-static) case (A = A 2 = Az> = 0) the expression above 
is identical to eq.(32) in j7|. 

Let us now assume that our system, initially at rest and in equilibrium, 
is perturbed and departs from equilibrium. Since uj is the fluid velocity 
as measured by our local Minkowski observer, then immediately after this 
departure we have 



uj m 



Co ^ 



and, from (|) and fllQl) , 



A « 



(43) 



(44) 



On the other hand, the following expression may be easily obtained for A 
from (|) and pi): 



A 
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[[P+ r) l-^2 + l-ufl 2 



UJ 



UJ' 



+ (p + P r ) UJ- 
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(45) 



or, evaluating this immediately after the departure from equilibrium, 

A = -87rre (l/+A)/2 [(p + P r ) uj] (46) 
Inserting (Q) and (H) into (^2|) , we finally obtain, using (|7|), (§) and @ 
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Note that boundary conditions 
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m s = M 







(47) 



(48) 



imply from ( 34]) and 



where M is the total mass (the mass parameter in the Scwarzschild metric), 

(m T )s = M + 



pc<j 



CO 



(49) 



Thus, although (m<r)£ = M in the static (and quasi-static) case and imme- 
diately after the departure from equilibrium, this is no longer true in general 
in the dynamic case. We shall discuss on eq.([|7]) in the next section. 



5 Discussion 

Let us asssume that our system consists of pure incoherent dust (P r = P± = 
0). Then the value of for a spherical region of radius r within E, im- 
mediately after the system departs from equilibrium, depends on three dif- 
ferent contributions. The first one, M (r/r^) 3 , is the gravitational mass of a 
static, homogeneous sphere, of radius r, within E. The second contribution, 
— 4/Tr 3 / r rE ^e^ +A ^ 2 /r 4 J (Jq f 3 p' 'df) dr , depends on the inhomogeneity of the 
energy-density distribution and will be positive if, as required by stability 
conditions, p' < 0. The last term 47rr 3 J r rs e A (p + P r ) Codr will contribute 
negatively in the case of collapse {Co < 0). Thus, density inhomogeneity 
tends to increase the Tolman mass within the sphere of radius r, if p' < , 
leading thereby (according to (|56])) to a faster collapse. This conclusion is 
also true in the case of slowly collapsing spheres, since ([47]) is identical to 
( P^D in |0 except for the last term containing Co. This last term, in turn, 



11 



tends to "stabilize" the system by decreasing (increasing) the Tolman mass 
in the process of collapse (expansion). 

In the case of an anisotropic fluid, the term Aur 3 J r rs e^ +A " 2 (2/r) (P r — P±)dr 
plays a similar role to the density inhomogeneity, with P r > P± (P± > P r ) 
corresponding to p' < (p' > 0). The fact that P r > P± (P± > P r ) leads to 



a stronger (weaker) collapse is, intuitively, in agreement with eq. (|16|) . This 
conclusion concerning the influence of local anisotropy is also valid in the 
slowly evolving regime 0. However, in this case, we only have information 
about the system at the moment it departs from equilibrium. A complete 
description of the evolution of the system, requires a full integration of the 
field equations. 

We have seen so far why and how the pace of collapse is affected by the 
density inhomogeneity and local anisotropy. Since the former is related to the 
formation of naked singularities, it could be speculated that local anisotropy 
might also lead to the formation of naked singularities. Of course, only a 
specific example could confirm (or discard) this suspicion. 
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